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Abstract 

Finite thermostats are studied in the context of nonequilibrium sta- 
tistical mechanics. Entropy production rate has been identified with the 
mechanical quantity expressed by the phase space contraction rate and 
the currents have been linked to its derivatives with respect to the pa- 
rameters measuring the forcing intensities. In some instances Green- 
Kubo formulae, hence Onsager reciprocity, have been related to the 
fluctuation theorem. However, mainly when dissipation takes place at 
the boundary (as in gases or liquids in contact with thermostats) , phase 
space contraction may be independent on some of the forcing param- 
eters or, even in absence of forcing, phase space contraction may not 
vanish: then the relation with the fluctuation theorem does not seem 
to apply. On the other hand phase space contraction can be altered by 
changing the metric on phase space: here this ambiguity is discussed 
and employed to show that the relation between the fluctuation theorem 
and Green-Kubo formulae can be extended and is, by far, more general. 



I am honored to be given this occasion to thank Jiirgen Frolich and Tom 
Spencer for the time they spent discussing with me and communicating their 
insights, projects and ideas. In particular I have drawn inspiration from the 
one among the two to whom I have been in closer contact, (JF), particularly 
at the time when he worked at IHES: his strong criticism of any "mathemat- 
ical deviations" from physical problems has been always extremely effective 
on my work. 
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1 Thermostats 



A mechanical interpretation of the entropy production rate in nonequilib- 
rium systems interacting with thermostats and possibly subject to external 
non conservative ( "stirring" ) forces has emerged from simulations and stud- 
ies on nonequilibrium statistical mechanics since the early 1980's, [H [2j [31 0] . 
It is interpreted as phase space contraction rate, as defined by the diver- 
gence of the equations of motion which we write symbolically x = f(x), i.e. 
= -Eid x J{x). 

General thermostats acting on a mechanical system, on which also ex- 
ternal non conservative forces may act, will be modeled as described in Fig.l 
and illustrated in the caption, [HE]: 




Xo, Xi, . . . , X n 



c 



(figl) 
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X 0i = -cWo(X ) - Ej cWj(X ,X,) + Ei(Xo) 
Xji = -diUj(Xj) - diUj(X ,Xj) - ajXji 



Fig.l: The 1+n boxes Co, T±, . . . ,T n contain No, Ni, . . . , N n particles, of mass m = 
1, whose positions and velocities are denoted Xo, Xi, . . . , X„, and Xo, Xi, . . . , X n 
respectively. The E denote external, non conservative, forces and the multipliers 
<Xj model the thermostats and are so defined that the kinetic energies Kj = ^X^ 
are exactly constants of motion with values Kj = ^NjkgTj, kg = Boltzmann's 
constant, j — l,...,n. The energies [To, Uj, Wo,j, j > 0, should be imagined as 
generated by pair potentials tpo, <Pj, fo.j short ranged, smooth, or with a singularity 
like a hard coreQ 

To imply Kj = in the above model the multiplier aj has to be aj = 
_ (Qj+Uj) w u eTe 

3Njk B Tj ' wnelt; 

Qj= f - Xj ■ d Xj W ,i(X , Xj) (1.1) 

is naturally interpreted as the heat ceded per unit time to the thermostat 
Cj. The phase space contraction rate, neglecting for simplicity 0{Nj l ), is 
computed from the equation in Fig.l to be 



1 Singularities of different type but care has to be excercised in formulating and by 
external potentials modeling the containers walls and for simplicity the assumption 
of smoothness (possibly in presence of a hard core) is made here. For the more 
general cases, like Lennard- Jones potentials, see [6]. 
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a(X) 



j 



Qi + Uj 

k B Tj 



(1.2) 



(each addend should be multiplied by the factor (1 — if 0(N- l ) is not 
neglected) . 

Of course o~(x) depends upon the metric used on phase space and on the 
density giving the volume element: both are arbitrary and Eq. (jl,2p yields the 
contraction rate for the Euclidean metric and density 1: i.e. for the Liouville 
volume. Because of such ambiguity o~(x) cannot have an immediate physical 
meaning. However its time average, and the fluctuations of its finite time 
averages over long time intervals, have an intrinsic meaning, independent 
of the choices of the metric and the density, [5] , at least if the motions are 
"chaotic" , see below. 

Some interesting concrete examples of the above systems are illustrated 
in Fig. 2 and Fig. 3. 




periodic boundary ("wire") 
mx = E — ax 



(fig.2) 



Fig2. A modern version of the classical Drude's model for electric conductivity. 

In Fig.2 a = ^r%t and this is an electric conduction model of N charged 
particles (N = 2 in the figure) in a constant electric field E and interacting 
with a lattice of obstacles (circles in the figure); it is "autotermostatted" 
(because Co and T\ coincide) in 2 dimensions. This is a model that appeared 
since the early days (Drude, 1899, [7]) in a slightly different form (i.e. in 
dimension 3 and with the thermostatting realized by replacing the -ax force 
with the prescription that after collision with an obstacle velocity is rescaled 

to |*| = V / 5^T. 

The thermostat forces are a model of the effect of the interactions be- 
tween the particle (electron) and a background lattice (phonons). This 
model is remarkable because it is the first nonequilibrium problem that has 
been treated with real mathematical attention and for which the analog of 
Ohm's law for electric conduction has been proved if N = 1, [8]. 

Another example is a model of thermal conduction, Fig. 3: 
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(fig-3) 



Fig3. A model for thermal and electric conduction. 

in which N§ hard disks interact by elastic collisions with each other and 
with other hard disks (JVi = N2 in number) in the containers labeled by 
their temperatures Ti, T^: the latter are subject to elastic collisions between 
themselves and with the disks in the central container Co; the separation 
reflect elastically the particles when their centers reach them, thus allow- 
ing interactions between the thermostats and the main container particles. 
Interactions with the thermostats take place only near the separating walls. 

If one imagines that the upper and lower walls of the central container 
are identified (realizing a periodic boundary condition) and that a constant 
field of intensity E acts in the vertical direction then two forces conspire to 
keep it out of equilibrium, and the parameters F = (T2 — T\,E) characterize 
their strength: matter and heat currents flow. 

The case T\ = T2, E ^ has been studied in simulations to check that 
the thermostats are "efficient": i.e. that the simple interaction, via collisions 
taking place across the boundary, is sufficient to allow the systems to reach 
a stationary state, [9]. 

Thermostat models similar to the above have been considered in the 
literature, [3j \TU[ lllj . A fundamental problem with the model in Fig.l is 
that it is not clear which detailed assumptions have to be made on the 
interactions to insure that almost all initial conditions evolve staying in a 
bounded region in phase space so that they can be expected to determine a 
stationary state. This can be called the "thermostat efficiency problem" and 
it is, for nonequilibrium, the analogue of the Hamiltonian stability problem 
in equilibrium, [12J . The experiment in [9] encourages the idea that the 
assumptions could be very general and fairly simple. In [13J a model like 
the one in Fig. 3 was studied but the confinement difficulty was avoided by 
requiring that also the total kinetic energy Kq in the central container was 
constant thanks to an extra thermostatting force — aoXo with a properly 
chosen oq. 

The model in Fig. 3 without thermostatting forces to keep Kj, j > con- 
stant, hence with a purely Hamiltonian evolution, has been carefully studied 
in [14] which also gives the clearest account on the so called "transient fluc- 
tuation theorem" improving and extending its earlier formulation in |15j . 
and obtains implicitly also a transient version of the result on fluctuation 
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patterns, analogous to the one derived earlier for steady states in |16j . 

In [TJ] there is also a careful analysis of the model in Fig. 3 with the aim 
of obtaining results for stationary states: stationarity is made possible by 
taking the thermostats infinitely large stressing the (formidable) problems 
that one should encounter in attempting a rigorous proof. 

In this paper (and in all my preceding ones) I have chosen to consider 
only finite thermostats with empirical thermostat forces and studied a few 
problems by introducing a single assumption, the chaotic hypothesis. 

2 Chaos 

Microscopic motions are in all possible empirical senses "chaotic". The 
paradigm of chaotic motions are the hyperbolic transitive systems: these 
are smooth systems whose evolution can be intuitively described by saying 
that each phase space point moves being seen by the comoving neighboring 
points as a hyperbolic fixed point. 

Another intuitive way to look at such systems is to say that the phase 
space points can be coded into sequences £ = (£i)^L_ 00 of symbols, say the 
digits 0,1,2, ...,q < oo, in such a way that the dynamics becomes the 
trivial shift of the sequence £, and all sequences which satisfy M^^ i+1 = 1 
represent one phase space point, M being a "compatibility matrix" with 
elements Mij = 0,1 which is transitive (i.e. Mfj > for some s). There 
may be ambiguities, i.e. different sequences may represent the same point, 
but this can happen on a zero volume set of points only, in close analogy 
with the familiar ambiguity in the representation of number by digits (where 
0.9999.. and 1.0000... are the same number). 

It is natural, at least for some |17} [TT1 |18j . to imagine that motions 
of complex systems, like gases or liquids, are chaotic in the simplest sense 
(which is also the strongest) of being hyperbolic transitive on the attracting 
sets (also called Anosov systems). The chaotic hypothesis, proposed in [TTj . 
see also [12J, reflects this remark. 

Chaotic hypothesis Attracting sets for mechanical systems are smooth 
surfaces on which motion is smooth, hyperbolic and transitive. 
This is an hypothesis that has to be considered in the same sense as the 
ergodic hypothesis for equilibrium statistical mechanics, [19]. Hence it might 
be at first disturbing. 

However disturbing assumptions are common in the literature and, nev- 
ertheless, are often fruitful. I just mention the assumption of periodicity 
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with equal period ( "monocyclicity" ) of the motions of mechanical systems: 
it was employed in the derivation of the second law from the action prin- 
ciple in Boltzmann, [20]: this assumption was considered also by Clausius, 
Maxwell, Helmholtz and was the basis of the early works on the mechanical 
interpretation of the second law, |21l I22j . At the time there must have been 
objections to such a bold assumption and someone must have declared, as 
it was done a little later about its modification into the ergodic hypothesis 
(and as it is done today about the chaotic hypothesis), that it is "a strong 
assumption as the periodicity (or ergodicity) hypothesis raises the question 
of which systems of practical interest are "periodic" ("ergodic"), since al- 
most none of them is actually such", see |23j . Similar statements can be 
found in the literature, even in good papers. 

Chaotic systems (in the above sense) admit a statistics (called SRB 
statistics, [241 [25| I26j). i.e. a probability distribution fi on each attract- 
ing set which, by integration, gives the average values of the observables 
G(x) on trajectories whose initial data x are randomly chosen, near enough 
to an attracting set, with a distribution with some (arbitrary) density: 



(G) = lim - V G(S j x) = / G(y)n(dy), with probability 1 (2.1) 

where x — ► Sx is a discretized time evolution map, obtained by timing 
observations on the occurrence of some selected event. Or in the (unphysical, 
yet customary and interesting) case of observations in continuous time 



(G) = lim i / G(S t x)dt= [ G{y)n(dy), with probability 1 (2.2) 

T^oo T Jo J 

where x — ► Stx denotes the evolution of the initial data x via the equations 
of motion, [27] . 

If motion is chaotic (i.e. hyperbolic, regular, transitive) the finite time 
averages 

7 =<G>r = ( 2 - 3 ) 

r i=o 

satisfy a large deviations law, i.e. fluctuations off the average (G) as large 
as r itself are controlled by a function C{l) convex and analytic in a (finite) 
interval (71,72), maximal at (G). This means that the probability that 
7 G [a, b] satisfies 



6 



P r ( 7 G [a,b]) ~ e ^ax [a , b] c(7) ; Va,6e (71,72) (2.4) 

and the interval (71, 72) is non trivial if (G 2 ) - (G) 2 > 0, [2H EH [29] . If £(7) 
is quadratic at its maximum (i.e. at (G)) then this implies a central limit 
theorem for the fluctuations of y/r(G) T , but Eq. (|3.4p is a much stronger 
property. 

Remarks: (1) The hypothesis holds also in equilibrium; if the system 
admits a dense trajectory in phase space it implies the classical ergodic 
hypothesis. 

(2) If the observable G has nonzero SRB-average it is convenient to consider 
instead the observable ^7) because it is dimensionless, just as in the case of 

(G) = it is convenient to consider the dimensionless observable 



G 



{G ' 2) 

(3) If the dynamics is reversible, i.e. there is a smooth, isometric, map / 
of phase space such that I 2 = 1 and I St = S-tl or in the discrete case 
IS = S 1-1 /, then any time reversal odd observable G, with non zero average 
and nonzero dispersion (G 2 ) — (G) 2 > 0, is such that the interval of (71,72) 
of large deviations for is at least (—1,1) provided there is a dense orbit 
(which also implies existence of only one attracting set). 

(4) The systems in the thermostats model of SecQ] are all reversible with / 
being the ordinary time reversal, change in sign of velocity with positions 
unaltered, and the phase space contraction o~(x) is odd under time reversal, 
see Eq. (ll.2p . Therefore if 0+ = (a) > it follows that the observable 

3=0 ^ 

has domain of large deviations of the form {—g,g) and contains (—1,1). 

(5) Since by Eq. (jl.2p a differs from e(x) = J2j>o T^F by t ne time derivative 
of an observable, it follows that the finite or infinite time averages a and of e 
have, for large r, the same distribution. Therefore the same large deviations 
function controls the fluctuations of p' above and of 

1 e(S^x) 

P=-}^ ' a + - \ a )SRB = {£)SRB, (2.6 

3=0 ^ 

and it has been shown, [111 130] and in a mathematical form in [31], that 
under the chaotic hypothesis and reversibility of motions on the attracting 
set, the function £(p) has the symmetry property 
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C(-p) = CO) - po'+j for a11 p e (-P)P) 



(2.7) 



and p > 1. This is the fluctuation theorem of [TT] (it requires a proof and 
therefore it should not be confused with several identities, see for instance 
[32] . with which, for reasons that I fail to understand, it has been often 
identified). The interest of the theorem is that it is universal, model in- 
dependent yielding a parameter free relation which deals with a quantity 
which has the physical meaning of entropy production rate and therefore 
has an independent macroscopic definition and is accessible to experiments. 

(6) Eq. (|2.7p is closely related to the theorem in p3], from which it dif- 
fers only because it deals with finite thermostats assuming the (strong) 
chaotic hypothesis, rather than dealing with infinite thermostats and as- 
suming (strong) ergodicity properties. In spite of the latter work several 
paper have appeared in the literature trying to get rid of the chaotic hy- 
pothesis without adding much (if anything) to the lucid discussion in [14] 
about the necessity of suitable assumptions in order to allow extending a 
transient fluctuation relation (which is an identity, requiring no assumption, 
on the full phase space, [32]) to a stationary one (which deals with proper- 
ties that hold on a subset of zero probability with respect to the initial data 
sampling). 

(7) The fluctuation theorem has several extensions including a remarkable, 
parameter free relation that concerns the relative probability of patterns of 
evolution of an observable and their reversed patterns, [331 HZ1 133] , related 
to the Onsager-Machlup fluctuations theory, which keeps being rediscovered 
in various forms and variations in the literature. 

3 Onsager reciprocity 

Another consequence of the fluctuation theorem are the Onsager reciprocity 
and Green-Kubo formulae for the infinitesimal deviations from equilibrium, 
[13] : the latter can be independently derived (in a simpler way) from the 
chaotic hypothesis and time reversal symmetry assumed only at equilibrium, 
[8], as it will be shown in the concluding comments, or as discussed from a 
somewhat different viewpoint in [35j. 

Here the aim is to show that the Green-Kubo formulae, hence Onsager 's 
reciprocity, can be regarded as the version at zero forcing of the fluctuation 
theorem for stationary states. 

In the case in which T\ = T<i = ■ ■ ■ = T and E = the system is 
in thermal equilibrium and its state is characterized by a probability dis- 
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tribution fiQ which invariant under the time evolution x — * StX generated 
by the equations in Fig.l. Setting x = (Xo, Xo, Xi, Xi, . . . , X n , X n ) it is 
remarkable that the distribution can be explicitly found, [3], as 

Mdx) = conste -/3(^(Xo) + E J>0 (^(X,)W(X ,X J )) + Xo(Xo)) 

• (n i>0 S[Kj - f A^Xn^o dX, dXj) (3.1) 
where (3 = (neglecting 0{NJ 1 ) for simplicity). Calling the "unper- 
turbed" energy H (x) = K (X ) + E,>o^(Xj) + E i>0 ^i( x o, X.j! and 
SiKjix),^) = S(Kj(Xj) - iNjTj)), Eq.JSIO, written more compactly, is 

Ho(dx) = const e - pHo{x) J] 5(Kj(x),T) dx (3.2) 

i>o 

which is a distribution in an ensemble which, for the system in Co, is equiv- 
alent to the canonical one (for No, Lq — > cxd, Nq/Lq = const if Lq is the side 
of the container). 

We now want to compare the average values of various currents that are 
switched on when E, the external forces, become non zero and the temper- 
atures of the thermostats become different: E / and Tj = T + Sj. More 
precisely we look for the relations between infinitesimal forcing actions and 
the corresponding currents, i.e. the susceptibility coefficients. 

The currents are related to the average values of the derivatives of the 
entropy production with respect to the forces (material currents) or to the 
temperature inequalities (heat currents). However the arbitrariness inherent 
in the phase space contraction generates interesting questions: for instance 
in the model in Fig.l the phase space contraction with respect to the Liou- 
ville volume is independent of the external forces E, see Eqs. (jl.ip . (|1.2p . so 
that 3eo~ = 0, while it is obvious that the external forces generate material 
currents, being non conservative. 

On the other hand even in equilibrium a thermostatted system exchanges 
heat with the thermostats: hence there is a production of entropy which has 
a zero average but which is not zero and equal to J2j T^f- 

It is therefore interesting to see, first, why in equilibrium [i.e. when the 
thermostats have all the same temperature and no external forces act) the 
SRB -average of vanishes, |36| . This is the case because the latter 

quantity is the derivative of (3H (x). In fact the derivative 0Hq is (3 times 



2 The kinetic energy of the thermostats is an additive constant and therefore is 
not explicitly written. 



9 



the work done on the system by the forces — aXj which equals J2j>o k^T 5 • 
This means that a(x) — (3Hq{x) = and therefore 

E 7% = ^o(x) (3-3) 

and the r.h.s. is a time derivative, hence it has time average. 

When the system is out of equilibrium (i.e. Tj ^ T and E / 0) the 
heat currents flowing into the thermostats divided by the temperature are 
generated by the entropy production rate jk(x) = dr k cr(x), while the mate- 
rial currents through the system are defined by minus the derivatives with 
respect to the acting forces of the work per unit time that they do, given by 
the corresponding derivatives of Hq. Thus given arbitrarily (3 the quantity 

a(x) = a(x) -l3H (x) (3.4) 

generates all currents up to a proportionality factor (here f3 is arbitrary). It 
can be computed as 

a(x) = £ 9j+B. - /3E • X - p ^(Qj + Uj) (3.5) 
i>o KB± i j>o 

because, by the equations in Fig.l, Hq = E • Xo — J2j>o a j-^j an ^ therefore 
H = E • X + E j>Q {Qj + Uj), see Eq.flTI]). 

Hence, discarding the time derivatives terms involving the Uj (parame- 
ters independent), the currents (at infinitesimal forcing) can be generated 
by the function 

-oM = gQ^-^)-E.X 05 L (3.0) 

The generating function ctq is °dd under time reversal and vanishes at 
equilibrium Tj = Ti, E = if T is chosen T = Tj] its derivatives with respect 
to the forcing parameters Tj , generate the heat and material currents 
and, at the same, time uq{x) differs from the phase space contraction by a 
time derivative. 

Note that a is also the phase space contraction of the volume in phase 
space, provided the latter is measured by the distribution 

JL{dx) = const e -^ H ^ x) Y[ 5(Kj{x),Tj) dx (3.7) 

j>0 
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In [37] a reversible system (like the model in Fig.l), has been considered 
in which the generating function for the currents gq vanishes for vanishing 
"thermodynamic forces" F = (Ti — T, . . . , T n — T, E\, . . . , E„) = and sat- 
isfies the fluctuation relation or, better, its extension in [37, Eq.(14)], has 
been considered. 

And it has been shown, [37] . that the products of the currents, generated 
by the thermodynamic forces, times (3 = , and defined by 

jm = d Fr a(x) = d Fr a (x) (3.8) 

are such that their averages J m = (j m )sRB have susceptibilities L mp = 
Jv which satisfy 

m y F=0 

1 f°° 

L mp = - dt({j m (Sf)j p (-))sRB - (jm}sRB(jp)sRB)\ F =0 ( 3 - 9 ) 

£ J —oo 

If the parameter j3 is properly chosen as mentioned above, i.e. (3 = 
(and only if so chosen), <7o will vanish when F = 0. Since a and cto differ by 
a time derivative they can be interchangeably used in the theory of the SRB 
distribution and therefore a$ satisfies the fluctuation theorem (because a 
does) ; the assumptions in the derivation in [37] apply and therefore Eq. (|3.9p 
yields Onsager's reciprocity L mp = L pm , and Green-Kubo formula. 



4 Work and entropy theorems. Comments 

(1) This extends considerably the results in [371113] removing the restriction 
on the phase space contraction to be the generating function of the cur- 
rents. The key is that the phase space contraction is only defined up to a 
time derivative of an observable and the generating function of the currents 
coincides with the phase space contraction only if the observable is properly 
chosen. 

(2) it is worth stressing that the extension of the fluctuation theorem needed 
to derive from it Onsager reciprocity is an important one: in [33] it was 
further extended to show (conditional reversibility theorem) that there is a 
simple relation between the probability that an observable F(x), even or 
odd under time reversal (for simplicity), follows in a time interval — r, r a 
"pattern" F(Stx) = ip(t) or the "reversed pattern" F(Stx) = (f(—t) provided 
the entropy production rate is fixed, [33]. A statement that can be colorfully 
quoted as .. relative probabilities of patterns observed in a time interval of 
size t and in presence of an average entropy production p are the same as 
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those of the corresponding anti-patterns in presence of the opposite average 
entropy production rate, |34| p. 476], or also [34|, p. 476], or ... it "suffices" to 
change the sign of the entropy production to reverse the arrow of time, or 
also ... a waterfall will go up, as likely as we see it going down, in a world 
in which for some reason, or by the deed of a Daemon, the entropy creation 
rate has changed sign during a long enough time, [121 p. 288]. We can also 
say that the motion on an attractor is reversible, even in the presence of 
dissipation, once the dissipation is fixed. Again variations of this property 
keep being rediscovered, see for instance [38] . 

(3) In the case of systems in contact with a single thermostat but in a 
stationary nonequilibrium because of the action of external forces the above 
analysis has also interesting consequences. The phase space contraction 

can be written as o~(x) = J2j>o %"^r J > as m Eq. (jl.2p . or by adding to it 
a time derivative as W(x) = o~(x) + (3Hq{x) which in this case is simply 

a(x) = E ^^ X ° = ~£~r- Therefore the fluctuation theorem, as pointed out 
by Bonetto: see [12, Eq.(9.10.4)], yields the following "work theorem" 

(e-P WT )sRB = h w= f - [ T W(S t x)dt (4.1) 

T Jo 

in the sense that the logarithm of the l.h.s. divided by r tends to as r —* oo. 
More generally the identity up to a time derivative of a, a, J2j>o fc~~r" anc ^ 

a o = Z)j>o( fcflT ~~ PQj) ~ 0^ ' -^-0) see Eqs. (13.3p - (l3.6j) . implies that, in the 
same sense as in Eq. (|4.ip , the finite time average P of any of the latter four 
quantities, denoted a, over a time r will satisfy 

(e- PT )sRB = 1, P= f - f T a(S t x) (4.2) 

r Jo 

which can be called an "entropy theorem": not only remarkable because 
it involves quantities that can be measured in experiments, [3], but also 
because here (3 can be taken arbitrary, so that Eq. (|4.2|) is an infinite number 
of relations. Actually if p = P/(o'o)srb the large deviations of p satisfy the 
fluctuation theorem symmetry Eq. (|2.7p . Note however that all such relations 
are special cases of the theorem in |16j . 

(4) A further alternative method to derive the Green-Kubo relations is 
in [39]. It will be illustrated, for completeness, in the simple case of a 
system interacting with only one thermostat and subject to several non- 
conservative external forces that will be proportional to parameters E = 
(Ex, . . . , E q ). Under the chaotic hypothesis the SRB average of the currents 
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Jm — J 

VSRB(dx)j m (x), with jm(x) = dE m cr(x) in presence of thermody- 
namic forces E, can be computed as the limit J m = lim^oo Ha(j m (St x )) > ^ 
Sf 1 is the map such that x — > Sf'x solves the equations of motion in presence 
of forcing forces with parameters E, and no is the equilibrium distribution 
Eq-dEE]), [391 M- Therefore 

p+oo ^ p 

J m = lim fx (j m (Sfx)) = dt— no(dx)J m (Sfx) 
t^°° Jo dt J 

+°° d f fj-o(dx) E E 

dt T+ / mcE ^ o(dS t x)j m (S t x) (4.3) 
o dt J fj, {dbfX) 



-■DC' 



d f no(dS E t x) 
o dt J iJ.o{dx) 

but by the comment preceding Eq. (|3.7p (considered with Tj = T) 

d no(dS® t x) E 

~i — = o-(b_ f x) (4.4) 

so that the chain of equalities in Eq. (|4.3I) yields 



X' 



J m = dt J a(S® t x)j m (x)[io(dx) (4.5) 
And taking into account that cf(x) = 0, if E = 0, and j m (x) = dE m cf{x) 



Lpm = d Ep J m \E=o = dt(J d Ep cr{S E t x) d Em cr(x) / u (dx))| E=0 

1 poo p 

= g J ^ dt ( J dE v <j{S* t x) d Em <j{x) A*o(dx))| B = = L mp (4.6) 

by time reversal invariance of the equilibrium distribution no, which is the 
Green-Kubo formula. 
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